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Abstract We review the problem of stripe states in strongly correlated systems, and some 
of the theoretical, numerical and experimental methods used in the last few years 
to understand these states. We compare these states to more traditional charge- 
ordered states such as charge density waves (CDW) and phase-separated sys- 
tems. We focus on the origin of stripe states as an interplay between magnetic 
and kinetic energy, and argue that the stripe state is generated via a mechanism of 
kinetic energy release that can be described via strongly correlated models such 
as the t-J model. We also discuss phenomenological models of stripes, and their 
relevance for magnetism and for the pinning of stripes by the underlying lattice 
and by disordered impurities. Recent experimental evidence for the existence of 
stripe states in different cuprate systems is also reviewed. 
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1. Introduction 

The problem of strong electron correlations in transition-metal oxides and 
U and Ce intermetallics has been a subject of intense research in the last 20 
years. This interest has been driven mostly by puzzling experimental findings 
in materials such as organic conductors and superconductors, heavy-fermion 
alloys, and high-temperature superconductors. These systems are character- 
ized by large Coulomb interactions, low dimensionahty, strong lattice cou- 
pling, and competition between different phases: antiferromagnetism, ferro- 
magnetism, spin density waves (SDW), superconductivity, and charge density 
waves (CDW). The strong interplay between different order parameters is be- 
lieved to lead to charge and spin inhomogeneities, and to a myriad of energy 
and length scales that makes the problem very difficult to treat with the meth- 
ods and techniques used for the study of Fermi-liquid metals. Although the 
problem of the quantum critical behavior of metals in the proximity of an iso- 
lated zero-temperature phase transition has been subject of much study and 
heated debate [1], our understanding of the problem of electrons close to mul- 
tiple phase transitions is still in its infancy. Here we will review both the ex- 
perimental evidence for the existence of certain inhomogeneous states called 
"stripe" states and some of the current theoretical approaches used to under- 
stand their origin, nature, and importance in the context of magnetism and 
superconductivity [2]. 

The formation of static or dynamic spin-charge stripes in strongly correlated 
electronic systems has been corroborated recently by several experiments, es- 
pecially in manganites [3, 4], nickelates [5-8], and cuprates [9-34]. The exper- 
iments span a large variety of techniques, from scanning tunneling microscopy 
(STM) [4,9-14], neutron (and x-ray) scattering [5,6,15-25], nuclear magnetic 
(and quadrupole) resonance (NMR and NQR) [7, 26, 27], muon spin rotation 
(yuSR) [28, 29], optical and Raman spectroscopy [8, 30, 31], transport [32], 
angle-resolved photoemission (ARPES) [33], and ion channeling [34]. 

Charge and spin modulated states, such as CDWs, Wigner crystals, SDWs, 
antiferromagnetism, and ferrimagnetism are common occurrences in many 
transition-metal compounds. These systems are characterized by an order 
parameter (such as the charge and/or spin density) that is modulated with a 
well-defined wave vector Q. Because of the modulations and the coupling to 
the lattice, these states usually present lattice distortions which are easily ob- 
served in diffraction experiments such as neutron scattering. In this regard the 
stripe states discussed here are very much like CDW/SDW instabilities, except 
that in CDW systems the ordered state is driven by a Fermi-surface instability 
(usually generated by nesting and/or Van Hove singularities), and Coulomb ef- 
fects are secondary because of good screening. The clearest example of such 
Fermi-surface effects occurs in Cr alloys, where the system undergoes a phase 
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transition into a CDW/SDW state [35]. The CDW and SDW ti-ansitions occur 
at the same temperature, and the charge order has a period that is 1/2 that of 
the SDW. The main difference between the phase transition in Cr and the stripe 
states to be discussed here is that the charge order in stripe systems occurs at 
higher temperature than the spin order [16]. Thus, on reducing the temperature 
the onset of charge order occurs first and the spins simply follow. In a weak- 
coupling analysis of Fermi-surface instabilities, this type of transition is not 
possible because of reconstruction of the Fermi surface due to the appearance 
of long-range order. Thus there are different energy scales for the charge and 
spin order that characterizes the materials discussed here. 

We should stress that the CDW and/or SDW instabilities in metallic systems 
are not trivial, and although we understand the basic mechanisms which drive 
these instabilities [36] our knowledge of their origin and effects on the elec- 
tronic degrees of freedom is far from complete. Systems such as transition- 
metal dichalcogenides [37] have a high temperature CDW transition with a 
very anomalous metallic phase and show in addition the phenomenon of "stripe 
formation" [38]. The stripes in these CDW systems are understood, however, 
because the CDW order is incommensurate with the lattice and therefore phase 
fluctuations of the CDW order parameter are allowed energetically. Local 
CDW phase-slips give rise to a filamentary stripe phase which, in fact, has 
a Fermi-surface origin. 

To understand the origin of CDW stripes one may consider the complex 
order parameter A for a CDW with incommensurate ordering wave vector Q. 
The free energy of the problem may be expressed as [39] 

F = Fo[5R(A)] + j a!r|^-^[|Q-(V-iQ)A|2 + At|QxVA|2]|(l.l) 

where Fq [x] is a minimal polynomial of x that respects the symmetry of the 
lattice and renders the free energy bounded from below. For a triangular lattice, 
for instance, it can be written as [39] 

Fo[3fi(A)] = a(r, T)^{Af + 6(r, T)^{Af + c(r, ^)3f^(A)^ (1.2) 

where the coefficients of the expansion are smooth functions of temperature. In 
particular for the quadratic term a(r, T) = oo(r)(r — Tjcdw) where Tjcdw 
is the transition temperature of the incommensurate CDW (ICDW) state. In 
(1.1) m* and k are parameters specific to the material under consideration 
and the derivative terms are written such that the free energy of the CDW is 
minimal when the ordering wave vector lies in the correct direction and has the 
correct wave length. In the case of an ICDW the order parameter is obtained 
by minimizing (1.1) to give 

^iCDwir) = Ao,7e^Q ■^ (1.3) 
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where Aq,/ = ^y2ao(TicDW — T)/{3co) for T < Tjcdw and zero other- 
wise. The parameter c(r, T) = cq = constant. Here we have assumed that the 
parameters in (1.2) may be expanded in a form such as 6(r) = 60+^1 cxpjfKj • 
r} where Kj are the shortest reciprocal-lattice vectors characteristic of the lat- 
tice symmetry. 

In the commensurate CDW (CCDW) case the wave vector of the order pa- 
rameter "locks" with the lattice so that its modulation becomes a fraction of the 
lattice wave vector Ki. Then one would replace Q in (1.3) by Ki/g, where q 
is an integer, and Aq by a value Ao,c which must be calculated from the free 
energy and depends in general on various coefficients of Fq in (1.2). The tran- 
sition temperature Tccdw is usually smaller than Tjcdw so that the generical 
behavior of the system consists of two transitions, first into an incommensu- 
rate phase and then into a commensurate phase [39]. In many systems the 
ICDW-CCDW transition does not occur and the system remains incommensu- 
rate down to very low temperatures [37]. 

In order to study the problem of the commensurate-incommensurate transi- 
tion, and the topological defects which appear due to incommensurability, one 
must generalize (1.3) to include phase fluctuations. These may be incorporated 
by writing the order parameter in the form 

A(r) = Aoe*^'^^"+'^^"\ (1.4) 

where 9{r) is the angle variable which determines the commensurability of 
the system: for the ICDW e{r) = (Q - Ki/q) ■ r, while for a CCDW 9 = 0. 
Because we are considering the simplest problem of a single CDW wave vector 
the problem becomes effectively one-dimensional if the variables are redefined 
in a new, rotated, rescaled, reference frame defined by s = (z, y) = |Q — 
Ki/g|r. In this case it is obvious that ^(r) = 9{x), and by direct substitution 
of ( 1 . 1 ) we find that the dimensionless free energy per unit of length relative to 
the commensurate case becomes 

Sf = Jdx^{[dJ{x)-lf + g[l-cos{qe)]} (1.5) 

where g is the coupling constant of the system and depends on the parameters 
of (1.1). The free energy in (1.5) describes a sine-Gordon model where the 
cosine term favors the commensurate state (9{x) = 0) while the gradient term 
favors the incommensurate state {9{x) = x). Thus 9 is the order parameter and 
the discrete symmetry 9 ^ 9 + 27r/qis broken in the ordered phase. There is 
therefore a critical coupling value gc that separates these two phases. However, 
it is easy to see that there are other solutions which minimize the free energy. 
In fact, variation of (1.5) with respect to 9 yields 

(P9 

—^=gqsm{q9) (1.6) 
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which has a particular solution 

eK{x) = - arctan , (1.7) 

where the boundary conditions are 9{x = —oo) = d9{x = —oo)/dx = 0. 
Notice that (1.7) changes smoothly from ^ = at a; = — oo to ^ = 2TT/q 
when X S> l/y^- In the context of the sine-Gordon model this is called 
a topological soliton or kink, while in the CDW literature [39] it is called a 
discommensuration. In general, the solution of (1.6) is given by [40]: 

2 

6{x) = - a^:csm[r]sn{^/gx/k,k)] (1.8) 

where 77 = ±1 and sn(n, k) denotes the sine-amplitude, which is a Jacobian 
elUptic function of modulus k. The sine-amplitude is an odd function of its 
argument u and has period AK{k), where K{k) is the complete elliptic integral 
of the first kind. By substituting (1.8) in (1.5) and minimizing with respect to k 
one finds that k 2y/g. In the Umit g ^ one has also k ^ and sn(n, k) = 
sm{u), whence 6(x) x as expected. On the other hand, when g ^ 1 one 
finds A; 1 and sn(it, k) « tanh(«), in which case (1.7) is obtained. For a 
generic value of g the solution has the form of a staircase. The plateaus in 
the staircase are multiples of 27r/(/ and correspond to regions where the CDW 
is in phase with the lattice, while in the transitions between the plateaus the 
CDW is not locked, leading to discommensurations. For large values of g the 
discommensurations are rather narrow and we find stripe states. These states 
have been observed experimentally in CDW systems [38]. 

In contrast to these CDW stripes, the stripe systems which we will discuss 
here have their origin in Mott insulators with very large Coulomb energies, 
whereas typical CDW/SDW systems are very good metals in their normal 
phase (Cr is a shiny metal while La2Cu04 is opaque and grayish). It is exactly 
the "mottness" of these systems which complicates the theoretical understand- 
ing of their nature. If we take seriously the analogy between dichalcogenides 
and cuprates we could think of the stripes as phase-slips of an incommensurate 
order parameter associated with the Mott phase. The primary question con- 
cerns the order parameter of a Mott phase. Antiferromagnetic order usually 
occurs in a Mott insulator but is essentially a parasitic phase (systems of spin- 
less electrons at 1/2-filling with strong next-nearest neighbor repulsion, and 
frustrated magnetic systems, can be Mott insulators without showing any type 
of magnetic order [41]). Unfortunately, the order parameter which character- 
izes the Mott phase in a finite number d of spatial dimensions is not known. 
Within the dynamical mean-field theory (d 00) the order parameter of the 
Mott transition has been identified with a zero mode of an effective Anderson 
model [42], but generaUzations for the case of finite dimensions have not been 
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established. The search for the order parameter which characterizes "mottness" 
is one of the important problems of modem condensed matter physics. 

The aim of this work is to summarize the current literature on the stripe 
phase in high- Tc superconductors. Although there is a consensus for the ex- 
istence of stripes in manganites and nickelates, no agreement has yet been 
achieved concerning the superconducting cuprates. Despite this controversy, 
the presence of stripes is now firmly established in La2-2:Sra;Cu04 (LSCO) 
[16, 17, 30-33]. In addition, recent experiments suggest that they may also be 
present in YBa2Cu307_5 (YBCO) [19-26,34], as well as in Bi2Sr2CaCu208+5 
(BSCCO) [12, 13]. In this work we discuss some of the different experimen- 
tal techniques which prove or suggest the presence of stripes in cuprates, and 
we present some theoretical ideas on the existence and relevance of the stripe 
state. 

This review is structured as follows: in section II, a survey of the theoretical 
derivations of the stripe phase as a ground state of models which are appropri- 
ate for describing doped Mott insulators is presented. In section m we discuss 
the experimental observation of stripes, and in section IV the role of kinetic 
energy. In section V we consider some phenomenological models, which pro- 
vide a means to go beyond the question of existence of stripes and allows one 
to predict measurable quantities. Finally, in section VI, we present the most re- 
cent experimental results for YBCO and BSCO, as well as the open questions 
and topics of debates. In section VII we draw our conclusions. 

2. Origin of stripes 

One of the main problems in condensed matter theory since the discovery 
of high temperature superconductors in 1986 [43] is related with the possible 
dilute phases of Mott insulators [44]. These materials have a large charge trans- 
fer gap, so at half-filling are insulating two-dimensional (2D) antiferromagnets 
well described by the isotropic Heisenberg model [45]. These are trademarks 
of "mottness" [46] and led Anderson [47] to propose that cuprates may be well 
described by a Hubbard model with large intra-site repulsion U. Later studies 
showed that close to half-filling and infinite U the model maps into the t-J 
model [48] 

H = -tY. PiaC3,aP + JY.^i-^j^ (1-9) 

<i,j>,a <i,j> 

where t is the hopping energy and J k, At^ /U <C i, ?7 is the exchange interac- 
tion between neighboring electron spins, Sj = c| ^da^gCi^fj (ci^a is the electron 
annihilation operator at the site i with spin projection a =T, i, and a" with 
a = x,y,z is a. Pauli matrix). In Eq. (1.9) P is the projection operator onto 
states with only single site occupancy (double occupancy is forbidden). Eq. 
(1.9) reduces trivially to the Heisenberg model at half-filling and describes the 
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direct interplay between the two main driving forces in the system, magnetism 
(characterized by J) and kinetic energy (characterized by t). 

The existence of a stripe phase in cuprates was first suggested in the context 
of Hartree-Fock studies of the Hubbard model close to half-filhng and at zero 
temperature (T = 0) [49]. This calculation is essentially analogous to the one 
used to study CDW/SDW transitions in metallic materials. For U < t, vertical 
stripes (parallel to the x- or y-axis of the crystal) were shown to be lower in 
energy [49-51], whereas for large U > t diagonal stripes were found to be 
energetically more favorable [51, 52]. The crossover from vertical to diagonal 
stripes was calculated numerically to occur at U/t ~ 3.6 [51]. Such calcu- 
lations have recently been generalized to finite temperatures, and the phase 
diagram was derived as a function of T and doping Uh, (Fig. 1.1 [53]). One 
important feature of these mean-field calculations is that they predict the for- 
mation of charge-ordered domain walls at which the staggered magnetization 
changes phase by vr. The magnetic order parameter is therefore maximal not at 
(tt/o, 7r/a) as in an ordinary antiferromagnet, but at an incommensurate vector 
(tt/o lb (5, vr/a) where 5 = ir/i <^ I and i denotes the charge stripe spacing. 
This incommensurabihty is an important feature of the stripe problem because, 
as we will demonstrate below, it leads to a reduction in the kinetic energy of the 
holes. These calculations, however, always predict that the stripe states possess 
a gap. The simple reason for this effect is that the only way in Hartree-Fock 
to reduce the energy of the system is by opening a gap at the Fermi surface. 
Furthermore, Hartree-Fock calculations in strongly interacting systems are not 
quantitatively reliable because they are unable to take fluctuation effects into 
account, and therefore should be considered only as providing some quahta- 
tive insight into the ground-state properties. They, however, do provide a "high 
energy" guide (snapshot picture) of the possible phases of the problem, and in 
fact they have been fundamental for the interpretation of certain experiments 
such as neutron scattering. 

Recent theoretical efforts have focused not only on determining the ground- 
state properties of the Hubbard, but also of the t-J model. Numerical tech- 
niques such as density matrix renormalization group (DMRG) [54], quantum 
Monte Carlo [55], and exact diagonahzation [56], have been applied to the 
2D t-J model with differing degrees of success. The main problem is that a 
strongly interacting problem like the t-J model is subject to strong finite-size 
and boundary-condition effects which are difficult to control. 

2.1 Numerical Studies 

Early numerical calculations on the t-J model have shown that for physical 
values of J/ 1 and close to half-filling there is a tendency for phase separation 
[57]. This phase separation can be pictured in the limit of i = (classical 
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Figure 1.1. Phase diagram in the plane of temperature T and hole concentration x (rih = 2x 
for LSCO) obtained by Machida and Ichioka [53] from mean-field studies of the Hubbard 
model. In the figure, N indicates the normal phase and C denotes the commensurate antifer- 
romagnetic phase. In addition, two other incommensurate phases exist, with vertical (VIC) 
or diagonal (DIC) stripe order. The VIC phase is metallic, but the DIC is insulating at zero 
temperature. 

limit) as a lowering of the system energy by placing all the holes together in 
order to minimize the number of broken antiferromagnetic bonds. This simple 
picture leads to separation into two distinct phases: a commensurate, insulat- 
ing region and an insulating, hole-rich region. It naturally overestimates the 
importance of the magnetic energy relative to the kinetic energy, and there- 
fore can be correct only when J ^ t. For finite values of t, the hole wave 
function delocalizes and this picture breaks down. The main question is for 
which values of J/t a phase separation may arise. Emery, Kivelson and Lin 
[57] found that phase separation can occur for infinitesimal values of J/ 1 suf- 
ficiently close to half-filling. These results, however, have been questioned in 
the light of more recent numerical data. There is no doubt that the t-J model 
undergoes phase separation for J S> t as all numerical calculations indicate. 
Close to the physical region of J < t, the current evidence for phase separation 
is weak, and so the issue remains controversial. 

DMRG calculations in large clusters [54] indicate the presence of stripe cor- 
relations in the t-J model. These studies, however, have been criticized on the 
basis of the special role of boundary conditions. Recent work on the Ising t- 
J model indicate that stripe formation does occur in this system, independent 
of the boundary conditions [58]. It was shown via non-perturbative analytical 
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calculations that minimization of the hole kinetic energy is the driving force 
behind stripe formation. This result has been confirmed by a number of nu- 
merical calculations in the t-J model [59], as well as in the t-Jz model [60]. 
Another important conclusion from these studies is that in the stripe phase the 
superconducting correlations are extremely weak. In fact DMRG calculations 
show that the stripe state is a CDW/Luttinger-liquid state with vanishing den- 
sity of states at the chemical potential, and thus is naturally insulating [54]. 
The DMRG calculations then support the idea that stripe correlations compete 
with superconductivity instead of enhancing it. This is consistent with exper- 
imental findings in Nd-doped LSCO, where the superconducting transition is 
reduced when static stripe order sets in [16]. 

Although the DMRG calculations were originally performed in an t-J model, 
essentially the same physics is found in the t-Jz model. The reason for the sim- 
ilarity between the t-J (where the spins are dynamical) and the t-Jz (where the 
spins are static) may be understood on the basis of the fluctuation timescales 
for each component in the problem. In the t-J model the spins fluctuate with 
a rate ^ h/J, while the timescale for hole motion is Th ~ h/t. When 
J/t < 1 (the physical regime of the model) one has Tg > t^, that is, the holes 
move "faster" than spins. In this case a Bom-Oppenheimer approximation is 
reasonable, since the spins have slow dynamics, and the two problems become 
essentially identical [56]. The advantage of working with the t-Jz model is 
that many of its properties are significantly easier to study both numerically 
and analytically. 

The introduction of next-nearest-neighbor hopping t' favors mobile d-wave 
pairs of holes for t' > 0, and single-hole excitations (spin-polarons) for t' < 
[61]. At t' = the stripe state is very close in energy to the d-wave pair 
state, and thus a small change in the boundary conditions or inclusion of small 
perturbations in the Hamiltonian can easily favor one many -body state relative 
the other. The quasi-degeneracy of different many-body states is an important 
characteristic of strongly correlated systems. Moreover, in real materials other 
effects may also be responsible for the selection of the ground state, that is, for 
lifting of the quasi-degeneracy. Indeed, by including lattice anisotropics, which 
arise in the low-temperature tetragonal (LTT) phase of LSCO co-doped with 
Nd, the stripe state can be easily selected [62-64]: Hartree-Fock calculations 
of the Hubbard model have shown that a very small anisotropy (on the order 
of ten percent) in the hopping parameter t is already sufficient to stabilize the 
striped phase, independent of the boundary conditions (open or periodic) [63]. 
Monte Carlo studies of the t-J model have also confirmed these results [64]. 
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2.2 Stripes and phase separation 

The problem of the formation of inhomogeneous states in a system with 
phase separation can be easily understood from a classical point of view by 
studying the Ginzburg-Landau free energy functional. Let ijj be the order pa- 
rameter of a system described by a free energy, F, of the form 

= a{T,xm^ + + (1.10) 

where V may be complex for a superconductor, a{T, x), /3{x), and 7 > 
are functions of the temperature T and some parameter x (such as doping or 
pressure). In the theory of second-order phase transitions, the term is 
neglected close to the critical line because /? > in this region of the parameter 
space. Here, however, we assume that P may be negative, and therefore this 
term is required so that the free energy is bounded from below. The critical line 
in the (T, x) plane is given by Tc{x) (as shown in Fig. 1.2) and we assume the 
existence of a quantum critical point (QCP) at x = Xa (that is, Tc{xa) = 0). 
Close to the critical line we introduce the parameterization: 

a{T,x) = Qo [T/Te(x) - 1] , 

P{x) = Po[Tcix)/T,{xs)-l], (1.11) 

while 7 is approximately independent of x and T. Notice that with this choice 
the parameter f3{x) is positive for x > Xg, signaling that in this regime the 
transition is of second order. However, (3{x) vanishes at x = Xs and becomes 
negative for x < Xg, indicating that the nature of the phase transition changes 
at small x. In fact the point {xg, Tc{xs)) is a tricritical point. 

For X > Xg the term is irrelevant close to the phase transition, and the 
transition is of second order depending on wether T is greater or smaller than 
Tc{x). Minimizing the free energy with respect to the order parameter yields 

for T < Tc{x) and x > Xg. Notice, however, that at x < Xg the parameter 
P{x) vanishes, and one must include the term. In this case, the free energy 
has two minima (instead of one) at the critical line indicating that the system 
has two phases, one with ip = (normal) and another with if) = ipQ (ordered). 
Minimization of F with respect to the order parameter provides the condition 
for the phase transition 

/?'(^) = ya(r*,x)7, 



T*{x) = re(x) + 



16q;o7 



Tc{xg) 



(1.13) 
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Figure 1.2. Temperature-doping phase diagram for a system close to phase separation. The 
symbols are explained in the text. 

The solution of these equations gives two critical lines, Ti{x) and T2{x) on 
Fig. 1.2. These Unes terminate at x\ and X2, and for x\ < x < X2 there is a 
coexistence region with two phases (normal and ordered). 

Long-range interactions are readily introduced by modifying the IV'I^ term 
in the free energy to 



In this case of fully phase-separated states the cost in electrostatic energy is 
too high and phase separation is frustrated to a finite length scale, ^ps, that 
depends on the coefficient of the IV'I^ term. The formation of finite droplets 
with if) = ipQ and size £ps is therefore more favorable than the separation of 
the system into two homogeneous phases with = and ip = ipQ. Stripes can 
also be generated in this model if one adds terms which break the rotational 
symmetry 



depending on whether the interaction with the lattice may be represented in 
terms of L = (plus sign) or L = 2 (minus sign) angular momentum states (s- 
and d-wave, respectively). In the L = case a checkerboard state is favored, 
but even a small d-wave term generates stripes along the x- or y-directions. 




(1.14) 




(1.15) 
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Disorder can also frustrate the phase separation as one may show by adding 
a "random mass" term 

Fd = j dr m(r)|V'(r)|2 (1.16) 

to the free energy, where m(r) is a gaussian variable with average zero and 
variance u. Using the replica-technique with n replicas (n — > at the end of 
the calculation) and averaging over disorder gives 

n n „ 

F = ^ - n ^ y dr\Mr)f\Mr)f, d-H) 

a=l a,b=l 

where F^^^ is the free energy without disorder and with n fields V'n- Notice 
that in the replica-symmetric case (tpn = for ^) the disorder generates a 
term of the order — ul^/il^ which decreases the effective value of /?, therefore 
reducing the value of Tc{xs). In a renormalization-group (RG) sense this term 
is relevant, and if the disorder is sufficiently strong it will bring the tricritical 
point to zero temperature (that is, the tricritical point becomes a quantum crit- 
ical point (QCP) and will completely destroy the first-order phase transition). 
However, in the ordered phase the system may still possess a coexistence phase 
with different values of the order parameter Once again we would have a sit- 
uation where the system forms droplets of the paramagnetic phase (tp = 0) 
inside the ordered phase. The size of these droplets depends on the strength of 
the disorder and it is easy to make them adopt a stripe conformation by adding 
a term of the form (1.15) that breaks the rotational symmetry. 

Although the phenomenology of the problem is quite clear, what is not so 
evident is how to apply this theory to the cuprates. Emery et al. [65] proposed 
that a model similar to the one discussed here (the Blume-Emery-Griffiths 
model) may be applied to the cuprates if one defines a pseudo-spin Si which 
takes the values Si = +1 and 5^ = — 1 on regions corresponding to hole-rich 
and hole-poor, respectively, whereas iSj = indicates a local density equal to 
the average value. In this case V'(r) is the coarse-grained version of Si and the 
above discussion is applicable. Notice, once again, that this model completely 
disregards the kinetic energy of the problem and can be applied only in the 
situation where t = 0. It is therefore not at all surprising that stripes appear. 
The inclusion of itinerant degrees of freedom is not straightforward. One of the 
main effects of the presence of itinerant degrees of freedom is the generation 
of dissipation which can change the dynamical properties (and exponents) of 
the system. This problem has been the object of recent intensive study in the 
context of quantum phase transitions when the coupling between the magnetic 
order parameter and the electrons is weak, and the electrons may be treated as 
a Fermi liquid [66]. In this case the electronic system serves as a heat bath for 
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the relaxation of tiie magnetic order parameter. In Mott insulators, the mere 
existence of a Fermi surface and a Fermi-liquid state can be questioned and 
thus it is clear that the weak-coupling formalism cannot be applied to these 
systems. The charge degrees of freedom cannot be modeled purely as a heat 
bath because their feedback effect in the magnetic system is very strong. 

It is interesting to compare the two mechanisms for stripe formation quoted 
previously. In one mechanism (represented by the Hartree-Fock calculations), 
stripes are long period CDWs arising from Fermi-surface nesting in a weakly 
incommensurate system [49-52]. There are four features which arise from a 
Fermi-surface instability: 1) the transition is spin driven, i.e., there is a sin- 
gle transition temperature Tc below which the broken-symmetry solution of 
the Hartree-Fock equations is stable; 2) in the low-T phase there are gaps or 
pseudo-gaps on the Fermi surface; 3) the spacing between domain walls is 
equal toir/x, where x is the hole concentration; 4) the high-T phase should be 
a Fermi liquid. 

In the case of stripes arising from Coulomb-frustrated phase separation, the 
situation is quite different [67]: 1) the transition is charge driven, i.e., local spin 
order between the antiphase domain walls can develop only after the holes are 
expelled from the magnetic regions. Ginzburg-Landau considerations indicate 
either a first-order transition, in which spin and charge order arise simultane- 
ously, or a sequence of transitions in which first the charge order and then the 
spin order appears as T is lowered [68]; 2) the stripe spacing is not necessarily 
a simple function of x, and there is no reason to expect the Fermi energy to lie 
in a gap or pseudo-gap; 3) a high-T Fermi-Uquid phase is not a prerequisite. 

3. Experimental detection of stripes in high-Tc cuprates: 
LSCO 

Although the first predictions for stripe formation in doped Mott insula- 
tors were made 13 years ago, not much attention was paid to these results 
in connection with high-Tc superconductors until 1995, when experimental 
data from neutron-scattering measurements in cuprates were interpreted con- 
sistently within a stripe picture [16]. Co-doping of cuprates has been extremely 
important for revealing the modulated charge states. However, the inclusion 
of co-dopant usually reduces Tc, raising doubts about the coexistence of su- 
perconductivity and the striped phase [16, 69]. The first experimental de- 
tection of stripes in the cuprates was achieved in a Nd co-doped compound 
La2-a;-2/NdySrj;Cu04. For y = 0.4 and x = 0.12, Tranquada et al. [16] found 
that the commensurate magnetic peak at Q = (vr/a, vr/a) shifts by a quantity 
5 = ir/£, giving rise to four incommensurate peaks. In addition, new Bragg 
peaks appear at the points {±26, 0) and (0, ±26), indicating that the charges 
form domain walls separated by a distance £, and that the staggered magneti- 
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zation undergoes a phase-shift of tt when crossing them. The position of the 
peaks indicates that the stripes are oriented along the vertical and horizontal 
directions, with a density of one hole per two Cu sites (quarter-filled). 

The reason why static stripes could be detected in this compound is based on 
a structural transition induced by the Nd atoms. Indeed, co-doping with many 
rare-earth species, including Nd and Eu, produces a buckling of the oxygen 
octahedra around the Cu sites and a corresponding transition from the low-T 
orthorhombic (LTO) to a low-T tetragonal (LTT) phase. The critical concen- 
tration of Nd needed to destroy superconductivity is a function of the charge- 
carrier density, i.e., the concentration x of Sr atoms. However, the buckling 
angle of the octahedra is a universal parameter: for tilts above a critical an- 
gle 9 ~ 3.6°, superconductivity is completely suppressed in these materials 
[69]. For the values of Nd co-doping y used in the first neutron-scattering ex- 
periments, superconductivity was still present, and the authors claimed that in 
their samples Tc ~ 5 K. However, coexistence of static stripe order and super- 
conductivity in LSCO is an issue that remains controversial, although recent 
experiments in La2Cu04+y have shown the coexistence of these two phases 
[70] in the same volume of the sample. 

Another important factor assisting the detection of charge stripes in LSCO 
systems by elastic neutron scattering was the selected doping concentration 
X ^ 1/8. The 1/8 anomaly was known since 1988, when electrical resistiv- 
ity measurements in La2-a;Baa;Cu04 (LBCO) were first performed [71]. A 
mysterious reduction of Tc was detected around x = 1/8, but the understand- 
ing of this phenomenon was possible only recently, in the light of the stripe 
picture. Indeed, Ba substitution also induces a structural transition, similar to 
rare-earth co-doping, which probably acts to pin the stripe structure, stabilizing 
the charge ordering, and hence reducing T^. Recently, Koike et al. have shown 
that the 1/8 phenomenon is common to all the cuprates [72] and that a similar 
effect must occur for x ~ 1/4 [73]. 

Though the first neutron scattering experiment was performed in a Nd doped 
sample with the "magic" hole concentration x = 1/8, further measurements 
on samples with x = 0.10 and x = 0.15 confirmed the existence of incommen- 
surate peaks in the spin and charge sectors, giving support to the stripe picture 
(Fig. 1.3 [16]). Moreover, systematic studies of superconducting LSCO sam- 
ples with a range of doping values x has been performed by inelastic neutron 
scattering [17]. The detected incommensurability is exactly the same as that 
obtained in co-doped samples (Fig. 1.4 [16]). Both elastic and inelastic neutron 
scattering, in addition to NMR [27], NQR, /iSR [28], Hall transport [32], and 
ARPES [33] measurements indicate that stripes are present in LSCO. A linear 
dependence of the incommensurability 5 as a function of the doping concen- 
tration X has been detected for a; < 1/8, indicating that the stripes behave as 
"incompressible" quantum fluids in this regime, that is, for 0.05 < x < 0.12 
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Figure 1.3. Experimental phase diagram obtained by Ichikawa et al. [16] for Nd-doped LSCO. 
Tch and Tm denote, respectively, the temperatures below which charge and spin ordering could 
be detected in this system by elastic neutron scattering measurements. The superconducting 
transition temperature Tc obtained by susceptibility measurements is also shown. In addi- 
tion, the structural transition lines from the low-temperature orthorhombic (LTO) to the low- 
temperature tetragonal (LTT) and to the low-temperature-less-orthorhombic (LTLO) phases are 
displayed. 

the hole density in each stripe is fixed (one hole per two Cu sites), and by in- 
creasing the amount of charge in the system one consequently increases the 
number of stripes and reduces their average separation £{x). Moreover, Ya- 
mada et al. [74] showed that in this regime Tc is also proportional to 6, i.e., 
Tc (X X (X 6 (X l/l{x). Above x = 1/8, however, the behavior of the system 
changes and 5 nearly saturates, indicating a transition to a more homogeneous 
phase. Recently, neutron scattering experiments were performed within the 
spin-glass regime, for 0.02 < x < 0.05 [18]. The result was surprising: the 
incommensurate peaks are rotated by 45° in reciprocal space, suggesting that 
the stripes are diagonal and half-filled, with one hole per Cu site, analogous to 
nickelate stripes [18]. However, this conclusion may be premature. Because 
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the incommensurate peaks are observed only in tiie spin, but not in the charge 
sector, other explanations of the phenomenon are plausible, such as the forma- 
tion of a spiral phase [75]. We will discuss this topic below in Sees. 4.2 and 5. 
A summary of available experimental data concerning the incommensurability 
is presented in Fig. 1.4. 

It is interesting to compare the two mechanisms proposed theoretically [49- 
53,67] for stripe formation in the light of the experimental results [16]. Charge 
order indeed appears before spin order, favoring the Emery-Kivelson proposal 
of frustrated phase separation [67], but the stripe separation clearly displays 
a linear dependence on the inverse of the hole density, as predicted by the 
Hartree-Fock analysis [49-53]. Concerning the stripe filling, the Hartree-Fock 
predictions are observed in the spin-glass regime, whereas the Emery-Kivelson 
proposition holds within the superconducting underdoped regime. However, 
recent slave-boson studies of the 3-band Hubbard model have shown that if 
the oxygen-oxygen hopping integral tpp is finite, quarter-filled stripes are more 
stable than half-filled ones [76]. DMRG studies of the t-J model found also 
that quarter-filled stripes are the lowest-energy configuration [54]. 
La2.xSrj(Cu0 4 




0.24 



Hole concentration (x) 



Figure 1.4. Summary of data concerning incommensurability 5 = 7r/£ as a function of doping 
concentration x. Data were obtained from neutron-scattering measurements by several groups: 
open and full small circles are from [18] and [17], respectively; dark squares are from [11 1], the 
grey one from [17], and the white one from [16]; large circles are from Refs. [17-22]. 

4. The single-hole problem: the role of kinetic energy 

The problem of a single hole in a 2D antiferromagnet has a long history 
and it is probably one of the best-studied cases of strongly correlated systems. 
This problem is by nature single particle because it deals with a single particle 
interacting with a complex magnetic environment. Physically this is realized 
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only in physical systems with extremely small carrier densities, and many of 
their properties can be related to the polaron problem [80]. Here we will not 
review the problem in any way (there are very good reviews on the subject) but 
we would like to stress the important role that the hole kinetic energy plays in 
determining the possible phases. Moreover, as we are going to show, the same 
physics may be responsible for the stripe phases in transition-metal oxides. 

The simplest limit of the t-J model with a single hole is the limit of J = 0, 
that is, U ^ oo. This limit was studied by Nagaoka [81], who showed that a 
single hole makes the system unstable toward a ferromagnetic phase. The ori- 
gin of ferromagnetism in this case lies in the minimization of the hole kinetic 
energy: because double occupancy is precluded, the kinetic energy conserves 
the spin projection, and there is no energy penalty for the formation of fer- 
romagnetic bonds, the kinetic energy is minimal when all the spins have the 
same direction. This process is essentially the same as that occuring in double- 
exchange systems such as manganites, where a ferromagnetic coupling be- 
tween the electron spin and a magnetic host produces a ferromagnetic state by 
minimizing the kinetic energy. For finite but small J {J/t <C 1) the same effect 
occurs, but instead of polarizing the entire plane of spins a single hole produces 
a ferromagnetic polarization cloud of size R: the hole gains a kinetic energy 
of order At — t/ {R/a)^ by being free to move in the ferromagnetic region but 
has to pay a magnetic energy cost of order J{R/a)'^^^ for the generation of 
a frustrated magnetic surface with ferromagnetic bonds. Minimization of the 
total energy of the hole indicates that the radius of the ferromagnetic region 
decreases according to R/a ~ (t/ ,7)^^^'^^^^ /{d — 1) as J/t increases (notice 
that for the case d = 1 this estimate always produces = oo for all values of 
t/J) [82]. 

As J/t increases the magnetic energy generated by the frustrated magnetic 
surface becomes too large, and R shrinks to zero. Thus, larger values of J/t 
lead to a change in the physics. A simple way to reduce the magnetic frus- 
tration is to reduce the frustrated surface of the spin configuration. Instead of 
a frustrated surface of misaligned spins it becomes energetically favorable to 
create strings of ferromagnetic pairs of spins due to the retraceable motion of 
the holes. It is clear that in this case the energy of the string grows hnearly 
with its size L, and therefore that the energy required to generate a string is ap- 
proximately JL. On the other hand, the hole kinetic energy changes from — t 
to a quantity of order —t + 1/ {L/af', and the problem is essentially equivalent 
to that of a single particle in a hnearly confining potential. The solution of this 
quantum mechanical problem is straightforward and minimization of the total 
energy shows that the size of the strings varies according to L ~ (t/ J)^/'^ inde- 
pendent of the dimensionality. Thus, on increasing J/ 1 the single hole case ex- 
hibits a crossover from a ferromagnetic polaron to the so-called spin-polaron. 
The confinement described here is not completely correct because the hole has 
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been considered as a semi-classical entity whereas in fact it is a fully quantum- 
mechanical object, which could undergo quantum tunneling over classically 
forbidden regions. This tunneling gives rise to "Trugman loops" where the 
hole can move diagonally but with a very small tunnel splitting (that is, very 
large effective mass) [83]. 

In any case, the true situation lies between these extremes and ferromagnetic 
polarization is concomitant with string processes. It is clear that the problem of 
the doped antiferromagnet centers on the compensation of magnetic frustration 
by reduction of the kinetic energy. Furthermore, the string mechanism provides 
a clear way to release kinetic energy, namely the retraceable motion of the hole. 
By incorporating both the kinetic energy (creation of ferromagnetic bonds) 
and the magnetic energy (generation of strings) one may understand how holes 
can move in a system with strong antiferromagnetic correlations. We have, 
however, discussed only the case of a single hole but for superconductivity it is 
important to understand the situation when the density of holes increases. The 
first crossover in these systems occurs when a finite linear density of holes 
(say N/L is finite but N/L"^ zero) is reached. The second crossover occurs 
when N/ becomes finite. Thus, in such strongly correlated systems one has 
always at least 2 crossovers: from single particle to ID and from ID to 2D. In 
the next section we will discuss the first crossover and show that it is related to 
the formation of stripes. 

4.1 Crossover from single particle to ID: stripes in the t-J 
model 

Consider an infinite antiphase domain wall oriented along one of the crys- 
tal axes directions of the system (Fig. 1.5). The cost in energy per hole to 
create such a state is J/2(n~^ — 1), where n is the linear density of holes 
along the stripe. However, the hole wave function is translationally invari- 
ant along the stripe and the kinetic energy gain due to longitudinal hopping is 
— 2isin(7rn)/(7rra). For J/t = 0.4 one may show that the energy is minimized 
for n = 0.32 with an energy Ej, k, — 1.255t which is larger than the energy of 
a hole in the bulk (spin polaron), given by Esp = — 2\/3t k, — 2.37t [56]. Here 
we have not included the transverse motion of the hole perpendicular to the 
stripe, which further reduces the kinetic energy of the system but gives also a 
finite width to the hole wave function. Using a retraceable-path approximation 
(but ignoring hole-hole interactions) one may calculate analytically the Green 
function for the holes [58]. 

Holes are confined to an antiphase domain wall by the potential generated 
by strings of overturned spins (that is, there is a linearly growing potential 
transverse to the stripe direction). One may show also that in this configuration 
the hole is actually a holon, i.e., it carries charge but no spin and any motion 
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Figure 1.5. Antiphase domain wall with one hole. Thick lines represent broken bonds, while 
dashed lines mark the position of the topological defect. 



of the hole away from the stripe produces a spinon, a particle with a spin of 
1/2 but no charge, of energy J. In the bulk the hole carries both spin and 
charge and therefore is a spin-polaron. Spin-charge separation is thus local, 
but not macroscopic [58]. Because of this effect, Trugman loops [83], which 
are responsible for hole deconfinement in the absence of antiphase domain 
walls, are not effective because the motion of holes away from the wall always 
produces an excitation of finite energy. 

One finds that for J/t = 0.4 and n « 0.3 the energy of the stripe state is 
Eq —2.5t, and therefore lower than the energy of the spin polaron. Fur- 
thermore, the width of the hole wave function has a value on the order of 3 
- 4 lattice units [58], and therefore extends a considerable distance from the 
antiphase domain wall, in contrast to the "cartoon" picture where the stripe has 
a width of only one lattice spacing [16]. Thus, considering a stripe as a com- 
pletely ID object is somewhat misleading because each hole may make long 
incursions into the antiferromagnetic regions. Moreover, it is clear from these 
analytic calculations that it is the single-hole kinetic energy which is responsi- 
ble for the stabilization of the stripe state. These results have been confirmed 
numerically by DMRG and exact-diagonalization studies [59]. In previous 
works semi-phenomelogical field theoretical models were proposed to explain 
the formation of anti- and in-phase domain walls and stripes [84-86]. How- 
ever, from the studies on the t-Jz model it becomes clear that stripe formation 
is a short distance problem (that is, it involve high energy states) and cannot 
not be properly addressed with the use of field theories that can only describe 
the low-energy, long-wavelength physics. 
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We should stress that we are discussing the ground state, that is, the lowest- 
energy stationary state and therefore the concept of "stripe fluctuations" refers, 
in this context, to excitations which are separated from the ground state by an 
energy of order t{J/t)'^^^ because of confinement in the transverse direction 
[58]. These results are consistent not only with the DMRG results for the t-J 
model for small J /t [59] but also with those for the t-Jz model [60]. As a 
consequence the stripe is metallic, in contrast to Hartree-Fock results which 
always produce a gap [87]. As in a Luttinger liquid [88] one expects that hole- 
hole interactions drive the system toward a CDW phase, which would become 
insulating in the presence of any amount of disorder [89]. Moreover, interac- 
tions cause the density of states vanish at the chemical potential following a 
non universal power of the interaction parameter [88]. 

4.2 Stripes, magnetism and kinetic energy release 

At half-filling, cuprates are antiferromagnetic Mott insulators with the Cu 
atom carrying a spin of 1/2. La2Cu04 is one of the most striking examples of 
a layered antiferromagnetic Mott insulator. The Neel temperature, instead of 
taking a value on the order of the planar magnetic exchange J (Ri 1500K), is 
approximately 300K. This occurs because of the low dimensionahty of the sys- 
tem and indeed, an 0(3) invariant 2D Heisenberg model may order only at zero 
temperature due to the breaking of continuous symmetry (Mermin- Wagner the- 
orem) [90]. However, the small inter-planar coupling, J_[_ (~ 10~^ J), stabi- 
lizes antiferromagnetic order at finite temperature. Magnetism in these systems 
is evident essentially across the entire phase diagram, although long-range or- 
der is lost with only 2% doping by Sr. Since the Sr atoms are located out of the 
Cu02 planes and their effect is to introduce holes, the doped holes are essential 
for the destruction of long-range order. More importantly, it is the minimiza- 
tion of kinetic energy of the holes which is the dominant mechanism for the 
suppression of magnetism. In the following, we discuss several different ways 
to understand the importance of hole motion in these systems. 

A first indication for the importance of hole kinetic energy is given by mag- 
netic measurements for x < 0.02 which observe the recovery of the magne- 
tization when the system is cooled below the so-called freezing temperature, 
Tp{x) [28]. The staggered magnetization Ms{x, T) vanishes at the Neel tem- 
perature T = Tn{x) and is a smooth function for Tpix) < T < T]s[{x). 
However for T < Tp{x) the magnetization seems to recover to the full value 
expected at x = and T = 0. This effect can be ascribed to localization of the 
holes after which they affect the magnetization only locally. Exactly where the 
localization of holes occurs remains unresolved. However, soft X-ray (oxy- 
gen K-edge) absorption experiments indicate that the holes are probably in 
the oxygen sites [91]. NMR experiments appear to indicate that holes would 
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localize preferentially close to the charged Sr atoms for electrostatic reasons 
[92]. However, because the system is annealed as temperature is reduced, it 
is quite possible that the unscreened Coulomb interaction between holes plays 
an important role in the localization process. If this is the case, localized stripe 
patterns may form even at low doping, although disorder effects from Sr dop- 
ing are very strong in this region of the phase diagram and one would expect 
any stripe pattern to be random in the CUO2 planes [62]. Recent neutron- 
scattering experiments at low doping find diagonal incommensurate peaks in 
the magnetic sector [18]. However, as no charge peak has yet been observed 
in the spin-glass regime, these measurements may be interpreted within the 
stripe model but the question remains open: the antiferromagnetic peaks could 
also be interpreted as the formation of a spiral phase [75]. In order to resolve 
this question unequivocally neutron-scattering experiments in samples heav- 
ily doped with spin-zero impurities such as Zn are required. If, as expected 
for a spiral state [75], the slope of the incommensurability as a function of Sr 
concentration x changes by a factor (1 — 2z), where z denotes the Zn concen- 
tration, the stripe hypothesis would be excluded in the spin-glass regime. 

Independent of the pattern of locahzation, holes may be localized either on 
the O or the Cu sites. If hole localization occurs at the O sites, one would 
expect a spin-glass phase to be observed at temperatures below Tp{x) because 
a localized hole at an O site liberates one spin (configuration p^) which frus- 
trates the antiferromagnetic coupling between neighboring Cu atoms [93]. If, 
on the other hand, holes are locaUzed on the Cu sites, the magnetization of 
the system would be reduced by one quantum of spin. Presumably, because of 
the delocahzation of the hole wave function between O and Cu atoms the two 
effects can occur simultaneously [48]. The key feature of these experiments 
is that they indicate the importance of the hole motion for the destruction of 
long-range order 

The rapid suppression of magnetism with hole doping can be contrasted 
with the slow suppression of antiferromagnetic order when Cu is replaced by 
a non-magnetic atoms such as Zn or Mg [94]. In this case, long-range order 
seems to be lost only at 41% doping, that is, at exactly the classical percolation 
threshold for a 2D Heisenberg system. As has been shown in recent theoreti- 
cal and numerical studies of the diluted quantum Heisenberg model, magnetic 
order seems to disappear only close to the classical percolation threshold even 
in the quantum system [95]. These results have been investigated experimen- 
tally and in fact the quantum fluctuations introduced by the dopants apparently 
are not sufficient to produce a quantum critical point (QCP) below the value 
of classical percolation [94]. Although the comparison between the problem 
of hole doping and Zn doping is not at first obvious, from the chemical point 
of view Zn introduces a static hole in the Cu plane because it has the same 
valence, but also has one extra proton. The situation on doping by Zn is there- 
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fore similar to the problem of holes localized at Cu sites. Indeed the recovered 
magnetization in the hole-doped case, M(x, 0), is very close to the value which 
one would obtain on replacing Cu by a density x of Zn atoms. Thus, further 
confirmation is obtained that the hole kinetic energy is the driving force behind 
the suppression of antiferromagnetic order in these systems. 

5. Phenomenological models: stripes and 
antiferromagnetism 

In the previous sections we have argued that the kinetic energy of holes is 
fundamental for understanding the stripe phenomenon in cuprates. However, 
dimensional crossovers are very difficult to measure experimentally. Phase 
transitions, on the other hand are easy to observe, because they produce strong 
effects in the thermodynamic properties. Because Mott insulators are usually 
antiferromagnetic one may ask if such dimensional crossovers in the hole mo- 
tion affect the antiferromagnetic phase? In the previous section we provided 
evidence that the kinetic energy of the holes is responsible for the destruction 
of antiferromagnetism in these systems when the hole concentration is of or- 
der of 0.01-0.02. How is it possible that such small doping levels can destroy 
a robust antiferromagnetic phase, with a Neel temperature which is of order 
300K? 

A possible way to understand the effect of the hole motion is to consider the 

formation of infinitely long stripes. The first obvious effect is a breaking of the 
spatial rotational symmetry. In the antiferromagnetic phase the spin rotational 
symmetry is also broken, indicating that both symmetries must be broken in the 
ground state of a striped antiferromagnetic phase. As a consequence the Gold- 
stone modes associated with the broken symmetries must carry information 
about them. For an ordinary (non-striped) antiferromagnetic phase these are 
spin- wave modes characterized by an energy dispersion E(k) = Cgk, where 
Cs is the spin-wave velocity, and a spin-stiffness ps associated with the twist 
of the order parameter [45]. In a striped antiferromagnetic phase the Gold- 
stone modes remain spin waves, but because of the broken rotational sym- 
metry their dispersion is different if the mode propagates along the direction 
of broken symmetry or perpendicular to it, i.e., the energy dispersion is not 
circularly symmetric and E{k\\, k^) = ^c^/c^ + c?^k\, where || and _L refer 

respectively to the directions parallel and perpendicular to the stripes. At wave 
lengths longer than the stripe separation i and energies lower than the first 
spin-wave gap due to the folding of the Brillouin zone, this kind of dispersion 
is guaranteed by the nature of the broken symmetries. 

In an ordinary antiferromagnet the spin-wave velocity is simply related to 
the lattice spacing a and the exchange constant J by = SJa'^. Thus, a 
striped antiferromagnet may be modeled simply by assuming that the only ef- 
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feet of the stripes is to introduce anisotropy in the exchange constants. Let us 
consider the case of a spatially anisotropic Heisenberg model with exchange 
constants Jx and Jy in the x- and y-directions, respectively, in which case the 
spin-wave velocities in each direction are given by [96] 



Microscopically, one may regard the stripes as causing local modification of 
the exchange across an antiphase domain wall from J to a value J'{< J). This 
alteration of J leads to a macroscopic change in the values of the exchange 
constants in the same way that the introduction of impurities in a solid leads to 
an average change in the unit-cell volume. 

To relate J and J' to Jx and Jy is not a trivial task. One possibility would 
be to solve the linear spin-wave theory for the striped antiferromagnet and 
calculate the derivative of the spin- wave energy at the ordering vector Q. For 
stripes with a separation of Ns lattice sites, this procedure requires the solution 
of Ng coupled differential equations. Besides being computationally intensive, 
the solution would not provide significant insight into the origin of J' and 
would just exchange one phenomenological parameter by another. 

The simplest theory describing a striped antiferromagnet is the spatially 
anisotropic non-hnear a model 



where n is a unit vector field. The symbols have been chosen to suggest the 
continuum limit of an underlying effective integer-spin Heisenberg Hamilto- 
nian on a square lattice [97]. The underlying anisotropy parameter is the ratio 
of the two exchange constants, or of the two velocities. 



The value of a characterizes the theory, but its exact dependence on micro- 
scopic parameters is not easy to derive. 

We proceed by making use of well-established techniques to analyze the 
behavior of the field theory described by the action (1.19) to predict the phys- 
ical properties of the system of interest. It is useful to reexpress (1.19) sym- 
metrically by a dimensionless rescaUng of variables x' = [ay ^i'^xlv, y' = 

(a)^/^yA(A ~ 1/a is a momentum cut-off), r' = j2{Jx + Jy)^/J^JySaT /Ti. 



Cl = 2S^a^Jy{J^ + Jy) 

cl = 2S^a'^J^iJ^ + Jy). 



(1.18) 




(1.19) 



(1.20) 
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The effective action (1.19) becomes 

Seff = j2^)j^ dr' J dx' J dy' (d^fif , (1.21) 

where /j, denotes x',y', and r', 

go{a) = ^coA/pO = [2(1 + a)/^/^ {aA)/S (1.22) 

is the bare coupling constant, cq = [2{Jx + Jy)^/J^Jy\^^'^{aS) /% the spin 
wave velocity and 

P° = sJjJyS'' (1.23) 

the classical spin stiffness of the rescaled model. The original anisotropy is 
now contained in the limits of integration. 

Notice that while a depends on the ratio Jx/ Jy, the spin stiffness depends 
on the product JxJy Thus, given a microscopic model where Jx and Jy are 
expressed in terms of microscopic quantities the field theory is well defined. 
Unfortunately no calculations yet exist for the microscopic form of these quan- 
tities, and certain assumptions are required concerning their behavior. If the 
spatial rotational symmetry is broken at the macroscopic level, that is, one has 
infinitely long stripes in the y direction, a simple choice would he Jy = J and 
Jx = aJ, whence 

P° = V^Pi (1.24) 

where 

PI = JS^ (1.25) 

is the spin stiffness of the isotropic system. This choice is valid only when 
the system is composed of a mono-domain of stripes [98]. If the system is 
broken into domains, in which case the rotational symmetry is broken micro- 
but not macroscopically the choice (1.24) may not be the most appropriate. At 
sufficiently long wave lengths the system is essentially isottopic and therefore 
the spin stiffness is the same in all directions, [96] 

P^s = PI, (1-26) 

which is obtained by choosing Jx = \foLJ and Jy = J j \fa. We note that these 
choices are essentially arbitrary and are based on quahtative expectations con- 
cerning the nature of the correlations at very long wave lengths. One may 
show that the choice (1.26) is appropriate very close to the antiferromagnetic 
phase (x < 0.02) where the breaking of the system into domains is quite prob- 
ably because of disorder effects [96]. However, at larger doping (x = 1/8) 
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the choice (1.24) is more appropriate because long-range stripe order is ob- 
served [98]. With these two parameterizations one may analyse the problem 
and calculate physical quantities for comparison with experiments. One gen- 
eral consequence of the anisotropy introduced by the presence of stripes is a 
growth of quantum fluctuations because of the reduction of effective dimen- 
sionaUty. In fact, by using large N methods and RG calculations, one may 
demonstrate that the effective spin stiffness is given by 

Ps{a) = p^s{(^) 

where 

gc{a) = \/87r2^Q;/(l + a) {in + Vl + a) 

+ v^ln[(l + v^rT^)/V^}~^ (1.28) 

is a critical coupling constant. Notice that (1.27) is reduced from its classical 
value p^g {a) for fixed anisotropy a, and that it vanishes at some critical value 
Oc, whence gdac) = 50- Hence, as a function of the anisotropy, the model 
exhibits quantum critical point where the system undergoes a quantum phase 
transition from an ordered Neel state to a paramagnetic phase. The loss of 
antiferromagnetic order at x = 0.02 can thus be considered as a consequence 
of the enhancement of quantum fluctuations due to the presence of stripes. 

As explained previously, these considerations are valid when the holes move 
along the stripes and modify the exchange constant across the stripe. However, 
when localization occurs at low temperatures (as observed in the recovery of 
magnetization in NQR experiments) the stripes essentially cease to exist and 
the system undergoes a phase transition into a spin-glass phase. The simplest 
way to understand this phase is to consider hole localization at the oxygen sites 
with consequent liberation of one spin 1/2 [93]. This spin frustrates the anti- 
ferromagnetic order, because the superexchange interaction of the O spin with 
the neighboring Cu spins is antiferromagnetic. This problem may be treated 
by considering the O spin as a classical locahzed dipole moment [75]. 

Finally, the theoretical results concerning the existence of stripes and the 
appropriate model for describing them can be summarized as follows: the pro- 
posal that the doped t-J model undergoes a phase separation is supported by 
variational arguments [77], diagonalization on small clusters [77], and Green- 
function quantum Monte Carlo calculations [78]. On the other hand, several 
quantum Monte Carlo calculations [55], series expansions [79], exact diago- 
nalization [56], and DMRG calculations [54] yield results contradicting these 
claims and supporting the stripe picture. In order to gain more insight into 
the problem, we begin by considering the simplest possible case, namely, the 
single-hole problem. 



5o(l) 
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6. Phenomenological models: Transverse fluctuations and 
pinning of stripes 

While much theoretical effort has been concentrated on determining whether 
stripes are the ground state of models such as the Hubbard and the t-J models, 
which are supposed to describe high-Tc superconductors, a parallel research 
direction has also developed which consists of phenomenological studies of 
the striped phase. In this case, one assumes the existence of stripes and dis- 
cusses further aspects such as their effect on the antiferromagnetic state (see 
previous section). Motivated by issues such as the static or fluctuating nature 
of stripes and the mechanism of stripe pinning, a phenomenological theory for 
the pinning of stripes has been developed. Zaanen et al. have related trans- 
verse stripe fluctuations to the restricted solid on solid model (RSOS) which 
describes the growth of surfaces [99]. In a simplified form of the model, the 
transverse kink excitations of stripes are mapped to a quantum spin-1 chain 
model [99], whose Hamiltonian is 



(1.29) 



Here, t has the role of a hopping parameter for transverse kinks, J controls the 
density of kinks and the D term represents a nearest-neighbor interaction of 
kinks. The spins take the values = 0, ±1, where +1 and —1 are associated 
respectively with stripe kinks and anti-kinks and = describes unperturbed 
(flat) segments. The full phase diagram for this problem was determined nu- 
merically by den Nijs and Rommelse [100] after earlier calculations by Schulz 
[101], who treated the spin-1 problem as two coupled spin-1/2 chains. We 
have recently reanalysed the calculations of Schulz [102] and derived the cor- 
rect phase diagram from this formalism, which agrees with the one obtained 
in Ref. [100] (see Fig. 1.6). Six different phases can be identified, depending 
on the values of the D and J parameters. If J is positive, the last term of Eq. 
(1.29) determines that 5*^ = and the stripe is straight (flat phase). If J is neg- 
ative, both values 5^ = ±1 are equally favorable with respect to the J term and 
the D term determines the value of S^. If D is positive, nearest-neighbor seg- 
ments prefer to be similar and therefore Sj^ and Sj^^^ will have the same sign. 
This gives rise to the ferromagnetic phase (diagonal stripes), with a sequence 
of kinks or anti-kinks. On the other hand, if D is negative, and 5^"*"^ prefer 
to have opposite signs, and the stripe will be bond-centered and flat, with a 
"zig-zag" shape (a kink follows an anti-kink and vice versa). In addition, both 
the flat and the bond-centered flat phases, which are gapped, can undergo a 
Kosterlitz-Thouless transition to gapless rough or bond-centered rough phases, 
respectively. The sixth phase, which was not identified by Schulz, corresponds 
to a gapped, disordered, flat phase (DOF), (Fig. 1.6). In contrast to the flat 
phase, this phase has a finite density of kinks and anti-kinks, which are posi- 



Charge Inhomogeneities in StronglyCorrelated Systems 



27 



tionally disordered, but have an antiferromagnetic order in the sense that a kink 
iSj; = 1 is on average followed by an anti-kink = — 1 (rather than another 
kink), with any number of = states in between them. The DOF phase is 
the valence-bond phase which is responsible for the Haldane gap. 




Figure 1.6. Sketch of the phase diagram for the spin-1 chain. The stripe configurations repre- 
sented by circles, and the corresponding values of are shown below. There are six different 
phases: 1) a gapped flat phase, corresponding to straight stripes {Sz = 0); 2) a gapless rough 
phase (spins equal to zero and ±1 distributed randomly); 3) a gapped bond-centered (BC) flat 
phase, which has a long-ranged zig-zag pattern (periodic alternation of Sz = 1 and Sz ~ —1); 
4) a gapless BC rough phase with a zig-zag pattern (antiferromagnetic correlations with dis- 
ordered 5z = ±1 but no Sz = states); 5) a diagonal stripe phase, corresponding to a fer- 
romagnetic state in the spin language; 6) a gapped disordered flat phase (DOF), where a kink 
{Sz = -1-1) is followed on average by an anti-kink {Sz = —1), but with some Sz = states in 
between. 

In the limit of negligible nearest-neighbor interaction (D ~ 0), the above 
model can be related to the t-J model [103] by considering the transverse dy- 
namics of a single vertical stripe in a frozen Neel background. The transverse 
dynamical properties of holes are described by the t-Jz model, and they move 
under the condition that the horizontal separation between two neighboring 
holes cannot be larger than one lattice constant. Notice that this condition does 
not restrict the motion of the stripe: it may still perform excursions very far 
from the initial straight-line configuration, but the line cannot be broken. In 
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this case the t-J Hamiltonian describing hole dynamics may be mapped onto 
a spin-1 chain Hamiltonian analogous to Eq. (1.29), but with D = Q [103]. 
A duality transformation of the initial quantum string Hamiltonian maps the 
problem onto one describing a ID array of Josephson junctions [104], which 
is known to exhibit an insulator-superconductor transition at {t/J)c = 2/7r^ 
[105]. This Kosterlitz-Thouless transition represents the unbinding of vortex- 
antivortex pairs in the equivalent XY model, which translates to a roughening 
transition for the stripe problem. In this way, the transition between the flat 
and rough phases along the vertical line {D = 0) of the phase diagram could 
be determined precisely by analytical means [104]. 

These models are related to the sine-Gordon model, and by studying the 
spectrum of the quantum string in a Hilbert-space sector of zero topological 
charge, the meaning of the transition in the "string" language may be clarified. 
At {t/ J)c the (insulating) pinned phase, which has an energy spectrum with a 
finite gap, turns into a (metallic) depinned phase where the spectrum becomes 
gapless [104]. This procedure allows the connection of two important and 
different classes of problems, namely the transverse dynamics of stripes in 
doped antiferromagnets and a system with the well-known properties of the 
sine-Gordon model. 

In all the models discussed hitherto, the pinning potential arises from the 
discrete nature of the lattice. However, the introduction of holes into the MO2 
planes (M = Cu or Ni) is not the only consequence of doping a Mott insulator 
such as La2M04. Doping an antiferromagnetic insulator also introduces disor- 
der into this material due to the presence of counterions, which act as attractive 
centers for holes. Doping with divalent atoms, such as Sr^+ produces quenched 
disorder, because the ionized dopants are located randomly between the CUO2 
planes. In contrast, doping with excess oxygen generates annealed disorder. 
Indeed, oxygen atoms have a low activation energy and remain mobile down 
to temperatures of 200-300 K. 

In order to account for the random pinning potential provided by the Sr 
atoms in nickelates and cuprates, one may add a disorder potential to the pre- 
vious phenomenological Hamiltonian (1.29) with D = Q. This allows a de- 
termination of the influence of both disorder and lattice effects on the striped 
phase of cuprates and nickelates. We consider the problem of a single stripe 
along the vertical direction confined in a box of size 21, where I denotes the 
sttipe spacing. The system is described by the phenomenological Hamiltonian 

n 

with t the hopping parameter, n„ the displacement of the n-th hole from the 
equiUbrium (vertical) configuration, p„ its conjugate transversal momentum, 
J the stripe stiffness, and Vn{un) an uncorrelated disorder potential satisfying 



-2tC0S r + J {Un+\ - Unf + Vn{Un) 



(1.30) 
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(T4,(tt)y„/(ti'))rf = d6{u — u')5ny, where denotes the Gaussian average 
over the disorder ensemble and d is the inverse of the impurity scattering time. 
Eq. (1.30) may be straightforwardly related to Eq. (1.29) by noting that S'f^ = 
Un+\ — Un and that the hopping terms 5^ and are connected to translation 
operators, which can be written in the momentum representation pn as r± = 
Q±ip„/h [104]. We are considering the lattice parameter a = 1. 

A dimensional estimate provides the dominant features of the phase dia- 
gram. At large values of the hopping constant t ^ J, the first term may be 
expanded as —2tcos{pn/K) ~ const. + t{pn/K)^. In the case of no impurity 
potential, Vn{un) = 0, hole dynamics is governed by the competition between 
the kinetic term t{knf', which favors freely mobile holes, and the elastic one, 
J{un+i — Un)^, which acts to keep them together. When the confinement 
is determined by the lattice pinning potential, the average hole displacement 
Un+i — Un is of Order 1 (the lattice constant is unity) so the wave vector kn ~ 1. 
A transition from the flat phase, with the stripe pinned by the underlying lattice, 
to a free phase is then expected at t/ J ~ 1. 

We now consider the opposite limit of strong pinning by impurities. In this 
case, the potential provided by the lattice is irrelevant and the typical hole 
displacement is on the order of the separation between stripes, 1/A;„ ~ Un+i — 
Un ~ By comparing the kinetic t{l/£)'^ and the elastic J(^)^ terms, we 
observe that a transition should occur sAt/ J ^ {l^. 

Indeed, by deriving the differential renormalization group (RG) equations to 
lowest nonvanishing order in the lattice and disorder parameters, one obtains 
a set of flow equations [89], which indicate that the transition from the flat 
(lattice-pinned) to the free phase occurs at (t/J)c = 4/7r^, and the transition 
from the disorder-pinned to the free phase occurs at {t/J)c = (36/7r^)£'^. The 
pinning phase diagram of the striped phase is shown in Fig. 1.7, in which 
5 = 1/21. 

By comparing these results with recent measurements on nickelates and 
cuprates one concludes that nickelates occupy the lower left comer of the phase 
diagram, i.e., they have static stripes which are pinned by the lattice and by 
the impurities. By contrast, cuprates are characterized by freely fluctuating 
stripes and so appear in the upper right corner. An appropriate treatment of the 
striped phase in cuprates must therefore include stripe-stripe interactions and 
the model becomes similar to that for a 2D fluctuating membrane [89, 106]. 
In this phenomenological framework, both nickelate and cuprate materials can 
be understood in a unified way, the difference between them being simply the 
parameter t/J which measures the strength of quantum fluctuations. 

Although the number of holes is intrinsically connected with the number of 
pinning centers, recent experimental developments show that it is possible to 
control these two parameters independently. Co-doping the superconducting 
cuprate material LSCO with Nd or Zn increases the disorder without modify- 
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Figure 1.7. Zero-temperature pinning phase diagram of tlie stripe piiase in tlie presence of 
lattice and impurity pinning. Three phases can be identified: a quantum membrane phase with 
freely fluctuating Gaussian stripes, a flat phase with the stripes pinned by the lattice, and a 
disorder-pinned phase [89]. 

ing the number of charge carriers [16, 72, 107]. On the other hand, growing the 
superconducting film on a ferroelectric substrate and using an electrostatic field 
as the control parameter allows the number of charge carriers in the plane to be 
altered for a fixed Sr concentration x [108]. This class of experiments consti- 
tutes an important step towards the control of a normal metal-superconducting 
transition. 

One may then investigate the stripe pinning produced by Zn and Nd co- 
doping [109]. The two dopants play fundamentally different roles in the pin- 
ning process. Nd, as with other rare-earth co-doping, induces a structural tran- 
sition which produces a correlated pinning potential trapping the stripes in a 
flat phase. The situation is analogous to the pinning of vortices by colum- 
nar defects or screw dislocations [110]. In this case transverse fluctuations are 
strongly suppressed, long-range order is achieved and thus the incommensurate 
peaks observed by neutron diffraction become sharper after the introduction of 
the co-dopant, as observed experimentally [16]. On the other hand, in-plane 
Zn- and Ni-doping provide randomly distributed point-Uke pinning centers, 
similar to oxygen vacancies in the vortex-creep problem. Within the model in 
which a stripe is regarded as a quantum elastic string, the effect of randomness 
is to "disorder" the string, promoting line-meandering, destroying the ID be- 
havior and broadening the incommensurate peaks. A perturbative treatment of 
the RG equations discussed previously show at the next higher order that this 
kind of pinning is relevant only in under-doped systems, in agreement with 
experiments [111]. 
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Finally, it is essential to go beyond the studies of transverse stripe exci- 
tations to consider the coupling between longitudinal and transverse modes. 
Longitudinal modes may be described by a Luttinger-liquid Hamiltonian, and 
the coupling between longitudinal and transverse fluctuations was investigated 
using bosonization [102]. One finds that a longitudinal CDW instabihty can 
arise if the stripe is quarter-filled and the underlying lattice potential has a zig- 
zag symmetry. This result has shed additional light on the connection between 
the formation of a LTT phase and the subsequent appearance of charge order 
in high-Tc cuprates (Fig. 1.3). 

Experimentally, the suppression of superconductivity in LSCO co-doped 
with Nd (Fig. 1.3), and also the upturn of the resistivity in the normal phase, are 
correlated with a structural transition from the LTO to the LTT phase. Indeed, 
at a; = 1/8 the charge ordering temperature Tco reaches its highest value and 
the superconducting temperature shows a local minimum. Neutron diffrac- 
tion experiments indicate that in the underdoped phase of LSCO the stripes are 
quarter-filled [16]. The formation of the LTT phase favors a zig-zag symmetry 
of the transverse stripe degrees of freedom. Thus, below Tltt the CDW insta- 
bility discussed above becomes relevant and stabilizes a bond centered string 
with zig-zag symmetry. If the stripe spacing is exactly commensurate, as at 
X = 1/8, a long-range-ordered CDW (Wigner crystal) can form, leading to the 
suppression of Tc- On the other hand, for incommensurate doping values, soli- 
tonic modes are present in the stripe which prevent long-range charge order 
Longitudinal charge order has hintherto not been observed. Nonetheless, the 
upturn of the in-plane resistivity below Titt suggests proximity to an insulat- 
ing phase inside the LTT phase of underdoped cuprates, which is likely to be 
the bond-centered zig-zag stripe [102]. 

7. Experimental detection of stripes in YBCO and 
BSCCO 

An important question at this point is whether charge stripes are pecuhar to 
the lanthanates or a generic feature of all the cuprates. The answer is not yet 
completely settled, and further experiments are required to achieve an unam- 
biguous conclusion. Measurements on YBCO and BSCCO compounds begin 
to provide a comprehensive picture. 

Inelastic neutron scattering experiments recently performed by Mook et al. 
on YBCOg.e [20], which corresponds to a hole doping of x = 0.10, detected 
a dynamical incommensurabiUty in the magnetic sector of 6 = 0.105 ± 0.01, 
which is exactly the value measured by Yamada et al. for the corresponding 
charge concentration in LSCO [17]. Later, magnetic incommensuration was 
also observed at 5 = 0.0625 in YBCOg 35, which corresponds to a doping 
X = 1/16 [21]. In addition, measurements of phonon broadening at a wave 
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vector consistent with the stripe picture (twice the magnetic wave vector) con- 
firmed the previous results for YBCOe.e and YBCO6.35 [21]. Eventually, a 
static charge order peak has been observed in a 21g crystal of YBCOg.ss, at 
a wave vector which is exactly double the dynamical magnetic incommen- 
suration, 25 = 0.127, as expected within the stripe picture [22]. Although 
the charge peaks are small, 6 orders of magnitude below the strongest crystal 
Bragg peak, their existence is undeniable. However, it should be emphasized 
that no charge order has been observed in YBCO6.5 and YBCOe.e, so the situ- 
ation is not yet resolved. 

Another important feature measured recently in YBCO is the ID nature of 
the stripes [23]. In a 4g crystal of detwinned YBCOe.e, one could observe not 
four, but only two incommensurate magnetic peaks (the second set of perpen- 
dicular peaks have nearly vanished, because the sample was almost completely 
detwinned). The results suggest stripes aligned along the 6— axis, in agreement 
with far-infrared spectroscopic measurements by Basov et al. [31], which in- 
dicate that the superfluid density is larger along the 6— direction. This behavior 
cannot be attributed to the chains, because in underdoped materials the missing 
atoms in the chains would inhibit the chain contribution to superconductivity. 

As a conclusion, one could state that neutron scattering experiments in 
YBCO at several different doping concentration [24] confirm the stripe pic- 
ture sketched for LSCO and reveal the universality of the previous results. 

Concerning BSCCO, the majority of experimental results are obtained from 
scanning tunnehng microscopy (STM) and spectroscopy (STS) [7-12]. The 
advantage of STM lies in its ability to measure simultaneously, with atomic 
resolution, both the surface topography and the local density of states (LDOS) 
of a material. The topographic image can be realized due to the exponential 
dependence of the tunneling current / on the separation of tip and sample. In 
addition, the differential conductance G = dl/dV, where V is the sample bias 
voltage, is proportional to the LDOS of the sample at the tip location. 

Low-T STS in BSCCO samples (T^ = 87K) revealed the existence of a 
large number of randomly distributed regions, with characteristic lengthscales 
of order 30 A, which have anomalous LDOS features. These features were 
initially referred to as quasiparticle scattering resonances (QPSRs) and were 
thought to be due to quasiparticle scattering from atomic-scale defects or im- 
purities [9], but later measurements performed in Zn-doped samples indicated 
that these inhomogeneities do not originate from impurities. The locations 
of Zn impurities were identified from the zero-bias conductance map and the 
LDOS map taken simultaneously at the same location showed no correlation 
between the intensity of the integrated LDOS and the location of the Zn impu- 
rities [10]. 

Spatial variations of the tunnehng spectrum and of the superconducting gap 
can be observed in pure as well as in impurity-doped BSCCO samples, and 
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seem to be intrinsic to the electronic structure. The gap ranges from 25meV to 
65meV and has a gaussian distribution [10]. The average gap is very similar to 
that reported previously from tunneling measurements. The spectra obtained 
at points with larger integrated LDOS exhibit higher differential conductance, 
smaller gap values and sharper coherence peaks, which are the characteristic 
features of spectra taken in samples with high oxygen doping concentration. 
This observation suggests the interpretation that the inhomogeneities may arise 
from differences in local oxygen concentrations [10]. 

Recent experiments in underdoped BSCCO indicate that these high-Tc ma- 
terials are granular superconductors, with microscopic superconducting grains 
separated by non-superconducting regions. By doping the material with Ni 
impurities, it was observed that the position of Ni atoms coincide with regions 
of small gap (A < 50 meV). In underdoped BSCO these small-gap regions, 
which have large G(A) are separated by percolative regions with large gap A 
and low G(A) [11]. 

Qualitatively new information is provided by STS measurements in mag- 
netic field up to 7.5T The quasiparticle states generated by vortices in over- 
doped BSCCO show a "checkerboard" pattern with a periodicity of four unity 
cells [12], in agreement with the charge periodicity expected within the stripe 
picture [16]. Indeed, the magnetic spatial periodicity previously detected by 
neutron scattering in overdoped LSCO in the presence [112] or absence [17] 
of a magnetic field is eight unit cells, exactly twice the charge spacing. Shortly 
after these measurements, the 4a periodicity was observed by Kapitulnik et al. 
in nearly optimally doped BSCCO without magnetic field [13]. Transforma- 
tion of the real-space data to reciprocal space showed a periodicity of 4a in 
the randomly distributed regions with anomalous LDOS, which was manifest 
in four distinct peaks in reciprocal space [13]. Two peaks corresponding to a 
periodicity of 8a along the diagonal due to the periodically missing line of Bi 
atoms in the BiO plane were also seen clearly, confirming the sample quality. 
However, later studies by Davis et al. have cast doubt on the stripe interpreta- 
tion: by analyzing the energy dependence of the wave vectors associated with 
the modulation, they have argued that the checkerboard LDOS modulation in 
BSCCO is an effect of quasiparticle interference, and not a signature of stripes 
[14]. Thus there is currently no agreement concerning the interpretation of the 
STM data. 

In summary, one may state that the presence of stripes is now firmly estab- 
lished in LSCO (static in Nd-doped LSCO, dynamical in pure LSCO). Recent 
neutron scattering experiments indicate that they are also present in YBCO. 
However, the magnetic YBCO stripes seem to be dynamical, whereas the 
charge ones are static. Concerning BSCCO, only STM measurements, which 
are just sensitive to charge and have the disadvantage of being susceptible to 
surface defects and pinning, were performed hintertho. They suggest that a 
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static charge stripe is present in this compound, although many further mea- 
surements are required for a definitive understanding. 

8. Conclusions 

Doped antiferromagnetic insulators have recently attracted a great deal of at- 
tention because many of the materials in this class exhibit novel and interesting 
behavior. The cuprates, for example, become metallic at low doping concentra- 
tion, and even superconducting at relatively high temperatures, whereas other 
systems, such as nickelates, show metallic behavior only at very high doping 
and are never superconducting. Manganites, on the other hand, can display 
the phenomenon of giant magneto resistance. Despite the different electric 
and magnetic properties exhibited by these compounds, a common perovskite 
structure connects them. Moreover, spontaneous symmetry breaking and stripe 
formation seems also to be a shared feature. 

Stripe-like ground states were first predicted from Hartree-Fock calculations 
of the Hubbard model [49-53]. Later studies of the t — J model also confirmed 
the initial results [54]. However, this theoretical work attracted significant at- 
tention only when Tranquada et al. detected the existence of static spin and 
charge order in Nd-doped LSCO [16]. The peaks measured by elastic neutron 
scattering were incommensurate with pure antiferromagnetic order, and sug- 
gested that the system had undergone a phase-separation into ID regions rich 
in holes (stripes), which were acting as domain walls in the staggered mag- 
netization. The fact that the magnetization changes phase by tt when cross- 
ing a domain wall has the consequence that the magnetic periodicity is twice 
the charge-stripe spacing. The associated incommensurability is therefore one 
half of the charge incommensurability, and the detection of both magnetic and 
charge incommensuration gave undeniable support to the stripe theory [16]. 
Later inelastic neutron scattering experiments in pure LSCO showed that there 
is a spin gap in these materials, but that incommensurate peaks can still be 
measured at rather low energies. The actual value of the spin gap depends 
on the doping. The incommensurabihty measured in these compounds is ex- 
actly the same as that obtained in Nd-doped LSCO, which is understandable 
considering replacement of La'^+ with Nd^"*" does not add charge carriers to the 
system, but induces a structural transition that helps to pin the stripes. The pres- 
ence of stripes in LSCO has been confirmed by several different experimental 
techniques. Inelastic neutron-scattering experiments [17-23] indicate that in- 
commensurate peaks are also detectable in YBCO, with a linear dependence 
of the incommensurabihty as a function of doping concentration similar to that 
observed in lanthanates. These results seem to indicate that the striped phase 
could be a generic feature of cuprates, instead of a peculiarity of LSCO. The 
requirement of very large samples imposed by the neutron scattering renders 
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experimental progress in the field very slow. Concerning BSCCO, at present 
only STM data is available and while these may support the existence of stripes 
[12, 13], the question remains open. It is important to note here that STM is 
a surface probe, whereas neutron scattering measures bulk properties. Thus, 
surface impurities could play a dominant role in STM experiments, but a sec- 
ondary one in neutron scattering. The years to come will hopefully show the 
truth behind all the controversies. 

From the theoretical side there has been steady progress in understanding 
stripe phases. Different approaches have been applied to the problem with 
differing degrees of success. Mean-field theories of different kinds [113-115], 
gauge theories [116], and quantum liquid crystal phenomenology [117] have 
been employed to describe the stripe state. In this review we focused primarily 
on studies based on the t-J model, where DMRG calculations have shown the 
presence of stripe phases [58, 60]. These studies have the advantage that rather 
few parameters determine the physics of stripe formation. Although this kind 
of approach describes very well the nature of the stripes, it seems to indicate 
that the stripe state is essentially insulating and extra degrees of freedom, such 
as phonons [118], may be required to explain the experimental data in these 
systems [119, 120]. 

Phenomenological models have also contributed to a significant evolution 
in understanding. Appropriate models for describing stripe fluctuations have 
been developed, and analogies with other known systems established [99, 103, 
104]. The effect of stripe pinning by impurities and by the underlying lattice, as 
well as the differing roles of rare-earth and planar impurities, has been clarified 
[89, 102, 104, 106, 109]. 

An important task remaining for the coming years is to show, both experi- 
mentally and theoretically whether and how stripes are connected to supercon- 
ductivity. Systematic investigation of the different ways to suppress supercon- 
ductivity may yield answers to this complex question. 
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